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Abstract 

We construct dijet differential cross sections at large momentum transfer, 
in which threshold logarithms have been summed to all orders in perturbation 
theory. This extends previous work on heavy quark production, by treating 
collinear singularities associated with hard, massless partons in the final state. 
The resummed corrections enable us to define, in the sense of factorization, the 
underlying color exchange mechanism. The influence of color exchange on the 
resummed cross section is contained in the eigenvalues and eigenvectors of an 
anomalous dimension matrix, which describes the factorization of coherent soft 
gluons from the hard scattering. The precise formulas depend on the partonic 
scattering angles and energies, as well as on the method used to define the 
jets in the final state. For cone dijets at fixed invariant mass, we find leading 
logarithmic corrections that, like those in the Drell-Yan process, are positive, 
and which grow with increasing dijet invariant mass. Other choices of dijet 
cross section can give, however, qualitatively different behavior, even at leading 
logarithm. 
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1 Introduction 



Factorized cross sections in perturbative QCD separate universal parton densities 
from process-specific factors that describe the perturbative hard scattering of the 
partons These factors, or hard-scattering functions, may include singular distri- 
butions, defined by their integrals with smooth functions, such as the parton densities 
relevant to the process. To be specific, for the production of a system of mass M, 
through the collision of partons of invariant mass squared s, the hard-scattering func- 
tion will contain, at nth order, terms as singular as (a"/n!)[ln 2ra_1 (l — z)/(l — z)]+, 
with z = M 2 /s. The limit z — > 1 is the edge of partonic phase space, or partonic 
threshold, at which no energy remains for QCD radiation. These singular contribu- 
tions often increase the cross section. The classic case is the production of Drell-Yan 
pairs, with M = Q the invariant mass of the electroweak vector boson produced by the 
annihilation of a quark pair. The possibility of controlling such singular distributions 
to all orders in perturbation theory was recognized early 0, and the full analysis, 
including nonleading logarithms was carried out |3|, |J. The formalism is closely re- 
lated to the even earlier resummation of logarithms in the transverse momenta of the 
pairs H]. Finally, it was noted that the resummation in logarithms of 1 — z can be 
performed, not only for the fully inclusive Drell-Yan cross section da/dQ 2 , but also 
for the cross section at fixed rapidity, da/dQ 2 dy |J. We will have occasion to recall 
and employ the explicit forms of these resummations below. 

In this paper, we shall derive threshold resummation for several physically-relevant 
dijet cross sections at large momentum transfer. We will find both similarities and 
differences compared to Drell-Yan and related cross sections, due to the color content 
of partons emerging from the hard scattering. One difference is that in jet production 
the hard scattering is itself a QCD subprocess, and involves color exchange. Another is 
that the kinematics at threshold depend upon the method used to define the event. We 
shall see that both of these features of jet production influence higher-order corrections 
in an essential way. The broad outlines of our reasoning, however, follow the derivation 
of resummation from factorization, as discussed in Ref . |7j . 

When a hard process is mediated by QCD, rather than color-singlet electroweak 
annihilation, the situation becomes, not surprisingly, more complex. Nevertheless, 
it is not difficult to recognize that the leading threshold logarithms associated with 
the initial state are the same for QCD processes as for Drell-Yan at each order in 
perturbation theory. This observation was made the basis for the first estimates of 
the effects of resummed perturbation theory in heavy quark production |8|, |10j- 
The extension to nonleading logarithmic corrections has been carried out recently 
for heavy quarks [1^, [13|]. At this level, intriguing effects begin to show up, in 
which the color structure of the underlying hard scattering influences the pattern 
of soft radiation near threshold. These results are related to the long-standing pro- 
gram of developing observable consequences of QCD coherence, in which properties 
of hadronic final states reflect color flow in partonic subprocesses |T3 . 



For dijet (as opposed to heavy quark) production, threshold singularities also differ 
from the Drell-Yan case because of the presence of the final-state jets, of momentum 
Pi, i — 1,2. When the cross section is evaluated at fixed dijet total invariant mass, 
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(pi +P2) 2 , these differences turn out to be at the level of next-to- leading logarithm 
in the singular distributions. As we shall show, however, when the cross section is 
evaluated at fixed (pi ■ P2), new leading-logarithmic contributions arise, which are 
negative. Depending on the partonic subprocess, these new contributions may even 
overcome the enhancements from initial state interactions, and produce an overall 
suppression of the cross section, relative to lowest order. 

The resummation of singular distributions at partonic threshold is important be- 
cause corrections very close to the edge of phase space may be numerically large, and 
taken at leading logarithm often (although not always) grow, rather than decrease, 



with the order []TT| . The resummed cross section also requires integrals over the ar- 
gument of the running coupling |3], |J, indicating enhanced sensitivity to soft gluon 
effects, and suggesting the necessity of incorporating nonperturbative corrections. In 
addition, the resummation of threshold singularities in QCD hard scattering requires 
a general approach to the interplay of color exchange at short distances with the 



pattern of gluon radiation into the final state [14]. We shall see that beyond leading 



logarithms it is in general not possible to separate initial from final state emission, 
and that, in fact, the color structure of the hard scattering directly influences the flow 
of energy into the final state, not only at low order, but to all orders in perturbation 
theory. 



As in Ref. [12|, we shall make strong use of the factorization properties of gauge 



theory cross sections at high energy. The basic complication in these arguments 
arises from the divergences of massless field theories in Minkowski space, in which the 
separation of long- and short-distance dynamics must be carried out relative to the 
light-like directions intrinsic to the problem. The principles and methods underlying 
the arguments below have been developed, for instance, in Refs. [[I], [T5|, |TB[ . We shall 
try, however, to make the discussion relatively self-contained, and to explain our use 
of technical results as they arise. 

We begin in Sec. 2 by identifying the cross sections, based on cone algorithms 
for jets, which we shall use as illustrative examples in this paper. In the following 
two sections, we go on to discuss the factorization properties of these cross sections 
(Sec. 3), and then to give interpretations for the various functions that appear in 
the factorized cross section in terms of field-theoretic matrix elements of nonlocal 
operators (Sec. 4). The pattern of these arguments is close to that used previously 
for resummation in heavy quark cross sections flT2"|, but now modified to treat the 



extra collinear divergences associated with the final-state jets. In Sec. 5, we organize 
the singular behavior of each of the component functions of the factorization formula. 
As in the case of heavy quarks, we identify an anomalous dimension matrix, which 
controls color-sensitive gluon radiation into the final stateQ. Combining these results, 
we give the resummed expressions for dijet cross sections in moment space. We 
observe that the form of our resummed cross section, Eq. (|7l|), does not depend on 
the details of the jet-identification algorithm, and applies to any infrared-finite dijet 
cross section at large momentum transfer. We conclude with a short summary, which 



1 In a forthcoming paper we will compute this matrix at one-loop order and discuss its diagonal- 
ization for the basic parton processes, including gluon-gluon scattering. 
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looks towards future work. 



2 Cone-based Dijet Cross Sections 
2.1 Definitions 

In this section, we introduce dijet cross sections, describing the inclusive hadronic 
production of a pair of jets, 

h A {p A ) + h B (p B ) — >Ji(pi) + J 2 (p 2 ) + X(k) , (1) 
at fixed rapidity interval, 

Av = 

2 \Pi p 2 

and total rapidity, 



A, = Im ^4 , (2 ) 



1 , pj+pj \ 

Vjj = - In — _ . 3 

2 \Pi + P2 J 

We consider the situation represented in Fig. |l|. We have two final-state jets, identified 
by cones of opening angles 8\ and 8 2 - The four- vector jet momenta Pi are defined as 
the total momenta flowing into the cones. In the limit of vanishing 8^, the p^ approach 
light-like momenta. 

To define the dijet cross sections, we must specify a large invariant, Mjj, which 
is held fixed. A natural choice is the dijet invariant mass, 

M]j = {p l+ p 2 )\ (4) 

but other specifications are possible. We shall see, in fact, that the nature of the 
resummed cross section depends critically on this choice. We shall illustrate this 
point with the alternate definition 

M 2 JJ = 2p 1 -p 2 , (5) 

the scalar product of the two jet momenta. In either case, large Mjj at fixed Ay 
implies a large momentum transfer in the partonic subprocess. 

As we integrate over allowable final states, the actual total momentum flowing 
into the cones may differ from a light-like vector by corrections proportional to the 
jet opening angle. To simplify our discussion, we shall assume that the cones are 
small enough so that all contributions proportional to 8i 1 may be neglected, but 
at the same time large enough so that 



a s (Q) In 



« 1 , (6) 



where Q is any of the hard scales of the cross section, typically the momentum transfer 

HE HI. 
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The introduction of cones removes all the final-state collinear singularities from 
the partonic cross section, which is then infrared safe, once the initial-state collinear 
singularities have been factored into universal parton distribution functions. 

We assume that the jets are identified by an algorithm that samples phase space for 
sets of particles flowing into cones of size oV In this search, the jet direction, and hence 
the center of the cone, is defined by the total momentum of the corresponding set of 
particles. The details of the jet identification algorithm are otherwise not important 
to our arguments. Nevertheless, it may be helpful to have a definite procedure in 
mind. 

To be specific, we define our cone jet algorithms as follows. Consider a final state 
1^0 = I {li}x)y consisting of the set X = {qi} x of particles with momenta q^ . For 
any subset of X, x C X, we define first a unit vector, 

n T = . (7) 



*l2j£x Qj 

Given h x , we define a new set, x'(n x ,5), consisting of those particles flowing into a 
cone of half-angle S, centered about n x , 

' - J Qk-n x \ 
x = \ qk '■ ^ > coso > . (8) 
I \Qk\ J 

In these terms, a jet is defined to be any set of particles, x C X, for which 

x = x', (9) 

i.e., for which the particles within the cone defined by the set constitute the set 
exactly. 

Of course, with this definition, a general final state may have many jets. To define 
our dijet cross sections, we take pi^ and to be the pair of jets of highest energy 
in the sample, consistent with some restrictions on their directions, implicit in our 
choices of ranges for yjj and Ay above. These algorithms are insensitive to emission 
of zero-momentum lines and/or rearrangements of momenta among collinear-moving 



particles. We may therefore expect them to be infrared safe [15, 1/]. Of course, as 



noted in [TjJ, the cancellation of virtual and real corrections fails when the final state 
reaches the edge of the region of phase space that defines a cone-jet cross section. 
(For a recent discussion of such momentum configurations in event shapes, see ||18|| .) 



This happens, for example, when a finite-energy line qi G x, reaches the boundary of 
the cone about n x , i. e., when 

= coso. (10) 

\Qi\ 

Such singularities, however, are restricted to a lower- dimensional surface in phase 
space, and thus do not spoil the finiteness of the inclusive jet cross section |T5| . 
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2.2 Factorization 



The standard factorized form of a hard dijet cross section is 
do-h A h B ^j 1 j 2 (S,5 1 ,5 2 ) 



dMjjdyjjdAy 



1 

52 



I dx A dx B (frf A /h A (xA,^)<Pf B /h B (xB,H 2 ) 

fA,fB=q,q,g' 



Pi ' P 

xH f A f B [ ^7^,a s (/i 2 )A,<5 2 



2 ' 



(11) 



where the explicit factor of 1/S 2 is introduced to make Hf A f B dimensionless. The 
incoming partons carry fractions Xa and Xb of the momenta of the incoming hadrons 
Ua and hs, respectively. These integrals are weighted by the nonperturbative but 
experimentally measurable parton distribution functions (densities), <pf/h- The fac- 
torization scale « separates the long-distance physics described by the parton dis- 
tributions from the short- distance hard scattering. The first argument, Pi ■ Pj/^ 2 , 
of the hard-scattering function, Hf A f B , represents all large invariants formed from 
the momenta of the incoming partons and final-state jets, Pi,Pj = xaPa,xbPb,Pi,P2, 
i 7^ 3- 

The perturbative factors Hf A f B are smooth functions only away from the edges 
of partonic phase space. The threshold for the partonic subprocess is conveniently 
parameterized by the variable z, 



xaxbS 



(12) 



where S = (j>a + Pb) 2 is the overall center of mass energy squared, and s = xa%bS is 
the corresponding quantity for the subprocess. At z — 1 (partonic threshold) there is 
just enough partonic energy to produce the observed final state, with no additional 
radiation. In general, H includes distributions with respect to 1 — z, with singularities 
at nth order in a s of the type 



ln m M 



1-z 



m < 2n — 1 



(13) 



defined, as usual, by their integrals with any smooth functions T{z) 



dz 



ln m (l - z) 
1 - z 



F{z) 



dz 



ln m (l - z) 



1-z 



\t{z) - m\ 



/ dz 
Jo 



ln m (l - z) 
1-z 



(14) 



All distributions of this sort have been resummed for Drell-Yan and more recently 
for heavy quark pair production cross sections at leading and nonleading logarithms 
0, [|, [T3. In the following, we extend this analysis to dijet cross sections. 

We can simplify Eq. (pTT|) , using the observation of 0, which applies to QCD 
as well as electroweak cross sections. According to Ref. 0, in computing the lead- 
ing power (1/(1 — z)) of H for z — > 1, we may treat the total rapidity yjj of the 



6 



produced heavy system (the dijets) as a constant, equal to its value at threshold, 
2/thresh = {l/2)ln(x A /x B ), taking -f A = Pa^+, Vb = Pb^— In tnis leading-power 
approximation, which includes all logarithmic corrections of the form of (p!3|), both 
leading and nonleading, we may write 

d(T hAhB ^j 1 j 2 (S,5 1 ,5 2 ) 1 ^ f 1 f , 2 \ 

dMjjdyjjdAy = P f £> J T ** J ) 

JJ yjJ y fA,fB=q,q,g 

, 2u/ M Jj\ J X i X A 

x 0f B /h B {xB^ )olz — I 5 [yjj - 2 ln ~J 

x Yl °fAfB^hh (l- z,— LL ,Ay,a s (fi 2 ),5 1 ,5 2 ) . (15) 

/l,/2 V ^ / 

In this form, we have introduced z, defined in Eq. (|T2p, as an integration variable, 
measuring the fraction of partonic energy squared going into the observed final state, 
in our case the jet pair. Its lower limit is 

z min = r = , (16) 

with Mjj defined as in Eq. (|j) or (|5|). With either choice (as noted above) z max = 1 
corresponds to partonic threshold. The simplified, dimensionless, hard-scattering 
function a, now depends only on 1 — z, Ay, the ratio of Mjj to /x, the coupling 
and the cone angles. Up to corrections of order Si, each jet evolves from one of two 
partons, fi and f 2 , emerging from the hard scattering, + — > /1 + /2, as indicated. 

To compute the hard-scattering function perturbatively, we turn to (infrared- 
regularized) parton-parton scattering, which obeys the same factorization properties 
as Eq. (ISj). The leading power as z — > 1 comes entirely from flavor diagonal distri- 
butions 4>f/f(x,fi 2 ) ||. Similarly, at leading power we may sum incoherently over the 
flavors of the partons fi, f 2 that fragment into the final state jets. Thus, for incoming 
flavors J'a, J'b, the partonic hard scattering functions may be written as a sum over 
partonic subprocesses /a + /# — > f\ + f 2 , as above, 

dMhdy.udAy = I ** J *Mu(**H ) 



dyjjdAy 

x Z ^/a/ b -/i/2 (1-2,— ,A«/,a,(/i 2 ),5i,5 2 ) , (17) 



}\j2=q,q,g 



with the same hard-scattering functions ^f A f B -*fxf 2 as above. To calculate <r to any 
order of perturbation theory, we construct the partonic cross section to that order, 
and factorize initial state collinear divergences into the light-cone distribution func- 
tions 0///, expanded in a s . The remainders give the perturbative expansion for the 
infrared-safe hard scattering function [|l|] . 
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Singular distributions of the sort (|T3| ) are now conveniently organized by taking 
a Mellin transform of the rapidity-integrated partonic cross section flTfl) with respect 
to r, Eq. (©, 



dr r 



N-l 



dyjj S 



2 do-fAfg^jij^S, 5 1} 5 2 ) 
dMjjdyjjdAy 



E h A /f A (N + l, /i 2 , e) 4> fa/fB (N + 1, /i 2 , e) 
f 

xcr f (N, Mjj/fi, a s (fi 2 ), 6 1} 6 2 ) , 



with a(N) = J* dz z N ~ 1 a(z), and 0(JV + 1) = ft dx x N 4>{x). To reduce clutter in the 
notation, we have denoted the set of 2 — > 2 partonic reactions /a + /b ~ > fi + f2, 
collectively by f . Divergent distributions in 1 — z produce powers of ln(iV), according 
to 



dz z 



N-l 



ln m (l - z) 
l-z 



-1 

m + 1 



In" 



+1 — + 0(\n m - 1 
N V 



N 



(19) 



In the following, we shall resum logarithms of N, from which the singular distributions 
of a may be reconstructed by inverting the transform. In the limit of large N, 
neglecting terms that decay as 1/N, we may replace iV + 1 by iV in the arguments of 
the distribution functions 6. 



3 Refactorization for Leading Regions 

We are now ready to construct a new, "refactorized" , expression for the partonic cross 
section, which generates all singular distributions of the form of Eq. (|T3| ) at partonic 
threshold. It will include, besides functions for the jet pair, a function that describes 
the soft, but still perturbative, radiation outside the cones, and which responds to 
the color flow in the hard scattering. The basic steps in the factorization process 
are quite similar to those already undertaken for Drell-Yan and heavy quark cross 
sections |T^[, but the kinematics of jet production is sufficiently subtle to warrant 
a description of how the general treatment must be modified to be applicable in this 
case. 



3.1 Leading regions and phase space 

We consider the purely partonic reaction in which flavors j a-, Jb, of momenta pa and 
Pb collide to produce the observed jets. The first step in writing a factorized form 
for this partonic cross section at threshold is to identify those regions of momentum 
space that contribute to the cross section at the (singular) leading power of 1/(1 — r) 
with r = Mjj/S = M 2 jj/(p A + pb) 2 - 

For a two-jet cross section these leading regions are illustrated in Fig. |2], which 
represents the cross section in terms of a cut diagram |l], 16|. Leading power 



contributions arise from momentum configurations in which on-shell lines of finite 
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energy fall into one of four "jet" sub diagrams. Two of these, labelled ipf A /f A and 
Tpf B /f B in the figure, represent the perturbative evolution of the incoming partons, 
and the motion of their fragments into the final state. Logarithms arise from regions 
in momentum space where internal lines of the diagrams approach the mass shell at 
the subspaces illustrated by the figure. In Fig. |2] and below, we shall work in a general 
axial (or temporal), n ■ A = 0, n 2 ^ 0, gauge. This simplifies the analysis somewhat, 
by insuring that collinear logarithms occur only internally to the jets fl], 15, y|. The 



soft subdiagram, S, generates only a single infrared logarithm per loop in this class 
of gauges. We emphasize, however, that our final results will be gauge invariant, and 
that it is possible to reexpress all of our arguments in covariant gauges. 

Returning to Fig. the ^'s give rise to a pair of active partons, whose scattering 
initiates the jet event. The two "short-distance functions" h and h*, to which the 
active partons connect, contain the effects of off-shell partons at the hard scatterings 
in the amplitude and its complex conjugate. The two remaining jets, J\ and J 2 , 
represent the fragmentation of the partons emerging from the hard scatterings into 
the observed jets. Finally, the function S represents the soft radiation, coupling 
incoming and outgoing hard partons. 

The leading regions pictured in Fig. |2| are identified by means of the analyticity 
and power counting techniques described in Refs. [15], 16]. Compared to the leading 



regions of the Drell-Yan process, which they generalize, they differ by the presence of 
the observed final-state jets, to which soft radiation may couple, and which require a 
sum over color structures for the hard scatterings. To anticipate, the total partonic 
cross section may be factorized into functions J«, h and S corresponding to these 
quanta. We shall have more to say below about how and why this may be done, 
and how to define the singular functions explicitly as vacuum expectation values of 
nonlocal, composite operators in QCD. 

We shall follow the procedure developed for the Drell-Yan process in Ref. Q by 
constructing the jet and soft functions in such a way that they absorb all singular 
dependence on 1 — r. This can be done in a straightforward manner, by matching 
the phase space of partons included in the ^'s, J's and S with the phase space of the 
underlying process. It is easiest to carry out this analysis in the center of mass (cm.) 
frame of the incoming partons. 

Consider first exact threshold for the two-jet process. At partonic threshold in 
the cm. frame, the jets J\ and J 2 carry equal and opposite-moving momenta, p\ and 
P2, 

pt+P , i = (Mjj^)=p t X+p%. (20) 

Now suppose we emit soft radiation of total momentum k outside the jet cones. The 
total squared invariant mass necessary for this radiation in the partonic process is 



^(l-r) 2 ), (21) 



S = ( PA + p B ) 2 = {p 1 +p 2 + kf = (pj + p 2 f + 2Mjjk + 0(M 
where we have used Eq. (^). In terms of the variable r, Eq. (fTB]), we then have, 

1-t = |^ + 0((l-r) 2 ) (Ml J = ( Pl+P2 f 
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1-t = ^f ! + ^ + ° ((1 - r)2) K = 2P!-P 2 )- (22) 

Thus near threshold, when Mjj = {pi + P2) 2 , the total invariant mass of the two 
jets, the variable 1 — r is given by the total cm. energy radiated outside the cones, 
up to corrections suppressed by a power of 1 — r, which we have agreed to neglect. 
On the other hand, for Mjj = 2p± -p2, the variable 1 — r receives contributions from 
the jet masses as well. This is a crucial difference, and leads to different behavior 
near threshold at the leading logarithmic level. The reason for this difference is that 
partonic threshold in the variable Mjj = 2p\ ■ p 2 is surprisingly more restrictive than 
in the variable Mjj = (pi +p 2 ) 2 - In the latter case, since we sum over all particles and 
momenta emitted into the cones, there is already considerable phase space available 
for the production of a pair of jets of the specified invariant mass (Mjj), total rapidity 
and scattering angles. The analogs of these thresholds in e + e~ annihilation are jet 
configurations in which all energy is emitted into two opposite-moving jet cones, with 
no soft radiation in the intervening directions. The sum over all such states has 
infrared, but not collinear divergences, as the latter cancel in the sum over allowed 
states. Then, when an infrared energy resolution SE is introduced, and soft radiation 
outside the cone is summed over up to SE, we find one logarithm of SE per loop [ |T6|| . 

For the case Mjj = 2pi ■ p 2 , Eq. ([22]) shows that r = 1 requires vanishing jet 
invariant mass in addition to energy flow only into the cones. Many fewer final 
states then contribute at threshold, specifically only those final states which consist 
of massless jets (at the cones' centers, according to the construction of Sec. 2.1). It is 
not surprising that large corrections arise as a result. Because they reflect a restriction 
in phase space, the additional corrections in this case are negative, corresponding to 
incomplete cancellation of collinear singularities. 

The refactorized cross section is pictured in Fig. |3], where the double lines repre- 
sent, as we shall see below, propagators associated with ordered exponentials of the 
gluon field in the jet directions. This decoupling of soft gluons from jet subprocesses 
is at the basis of proofs of factorization for single-particle inclusive cross sections ||TP| 



and of factorization for inclusive cross sections pO], |21J . We shall discuss the technical 
justification for the form of factorization proposed here in the next subsection. The 
physical basis of jet/soft factorization, however, is the inability of soft gluons to re- 
solve the internal substructure of the jets. The effect of all attachments of soft gluons 
to the jets is summarized simply by gauge rotations on the partons that connect the 
jets to the hard scatterings [[[], The gauge rotations may be represented by the 
ordered exponentials of the gauge fields that we shall encounter below. 

Given the structure in Fig. [| we may write the partonic cross section daf A f B ^j 1 j 2 
in a preliminary factorized form, as the product of functions which describe the short- 
and long-distance processes in the scattering, 

da }AfB ^ JlJ2 {SMM _ 1 mrf)( M » AyjQs(/ ^ 

(f) /Mjj (1 - t)Mjj 2 \ 

xKU , , Ay, a s (/r), <Si, S 2 



olMjjdAy 5 2 U 
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+0(1 -t) 



(23) 



Here H IL includes the hard scatterings. The overall rapidity, yjj, is integrated over 
(although for r near 1, it is always small in the cm. system). At the same time, 
both the long- and short-distance factors depend upon the rapidity difference Ay, 
or equivalently, the cm. scattering angle. The indices /, L describe possible color 
structures at the hard vertices, for the scattering amplitude and for the complex 
conjugate amplitude. We may further factorize the fully short-distance function into 
contributions from the amplitude and its complex conjugate, 

Hfl (^Ay,<f))=h*f {^,Ay,a s {^ h? {^,Ay,a s {^ . (24) 

As r — ► 1, only purely virtual loops can carry momenta of order Mjj. These are the 
contributions kept in Hn- The color structures L and /, of course, depend implicitly 
on the flavor content of the hard scattering, denoted by f. For example, when the 
underlying partonic hard-scattering is qq — > qq, there are only two possible color 
structures, which may be chosen as singlet and octet in the s or i channel. Similarly, 
the color-dependent factors h depend as well upon the spin of outgoing scattered 
lines. The dependence on spin, however, follows the same pattern as dependence on 
flavor, and we shall not exhibit it explicitly. We always assume unpolarized collisions. 

It is also useful to factorize the long-distance part of the cross section, Kli, sep- 
arating partons that are collinear to the incoming quarks from those that are in the 
final-state jets and those that are soft and "central" in rapidity. Referring to Eq. 
(f2~2|), the sum of the weights associated with the functions is (1 — r), so that the 
convolution is of the form 



K LI 



{{) (Mjj (l-r)Mjj 2 \ 

V At fJL J 

= \ dwA dwB dwi dw2 dws S(l — r — wa — wb — w\ — w 2 — Ws) 
Jo 

v / ( MjJ ( 1\ \ I ( MjJ ( 1\ \ 

x WUISa \ w a, ), e I ipf B /f B \ w B , , <W ), e I 

{i )fW S Mjj 2 \ 

x s li[ ,Ay,a s (fi ) J 



+0(l-r). (25) 

Dependence on the gauge vector is present in each factor, but has been suppressed. 
The five functions ipA, ipB, and S^j are all convoluted together in terms of 

the weights that each function contributes to the final state. This is possible because 
the weights identified in Eq. (|22|) above are additive in the particles of the final state, 
and because the factorization implies that there is no interference between final-state 
particles associated with the different functions. Let k$ be the energy of particles 
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emitted by the soft function S outside the cones, and similarly for the four other 
functions. Then the lack of interference implies that 



k° = k\ + k% + Jfe? + k Q 2 + k° s , (26) 

with k° the energy in Eq. (|22|). The complete contributions of each of the functions 
to the weights are given by 



W A = -rr~ W B 



Mjj " Mjj 

p 2 1 + 2Mjjk° 1 _ p 2 2 + 2Mjjk: 



w l — T7o w 2 



2 



Mh ^ Mjj 



= (27) 



M 



for Mjj = 2pi • p2, and by 



././ 



2/c^ 2/c^ 



w l = T7 — w 2 ~ 



Mjj Mjj 

= (28) 

for Mjj = (p\ +P2) 2 - Note that the product i^a^bJi^S of jet and soft functions 
behaves as (waWbWiW2Ws)~ 1 when 1 — r vanishes || 0, which is how logarithmic 
enhancements arise. We will be able to neglect contributions to pf except from J&i 
because we neglect corrections proportional to 5{. Finally, the dependence on the jet 
opening angles is included entirely in the jet functions, up to corrections proportional 
to the 5j. Although it is possible to do so, we shall not attempt a resummation in the 
opening angles (recall our assumption, Eq. @). Before exploring the consequences 
of these factorized expressions, let us discuss their justification. 



3.2 Relation to factorization for inclusive processes 

The details of the proof of the factorization of soft quanta from a jet depend on 
whether the jet in question is "initial-state", as the ipi in Fig. [3], or "final-state", 
as the Jj. Factorization is somewhat simpler in the latter case, and was discussed 



in Ref. |B| in the context of single-particle inclusive annihilation, and in Ref. in 
the context of transverse momentum distributions. The arguments are essentially 
identical in this case. The factorization of soft interactions from initial-state jets 
is an essential ingredient in proofs of factorization for inclusive processes, and was 
extensively discussed in Refs. [ P0| , ETJ in this connection. 

In factorization proofs for inclusive cross sections, the goal is to show that soft 
interactions cancel, but the factorization of soft quanta from jets is an important 



subsidiary result in this demonstration |^TJ. It is necessary, however, to distinguish 
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the factorization of soft quanta from partonic jets in Eq. Qzq), in which z is fixed 
near unity (near threshold), from their factorization in the inclusive case, where z is 



freely integrated over. Let us recall a few of the relevant arguments of Ref. [^, 
as they apply here. First, the difference between initial-state and final-state jets lies 
in the singularity structures of the jets, in terms of the momentum components of 
soft quanta. For final-state jets, all poles are in the same half-plane, corresponding 
to final-state interactions only. This simplifies the analysis considerably. Indeed, a 
short argument based on contour deformation allows us to use the eikonal, or "soft" , 
approximation for soft gluon momenta within jets, and to replace all the details of 
the jets' interactions with soft gluons by ordered exponentials [| . 

For initial-state jets, however, contour deformation is severely restricted, on a 
graph-by-graph basis, as poles from final-state and initial-state interactions conspire 
to pinch the momentum integration contours in regions for which the eikonal approx- 
imation fails [IT!|. It is, in general, only after the sum over final states that final-state 



interactions cancel, the pinches disappear, and soft quanta may be factored. When z 
is fixed near unity, or when moments are taken, final states are not all summed over 
with the same weight, and we may not, in general, factor soft quanta from jets in the 
same manner as in Ref. ||21| , without leaving over apparently nonfactoring remain- 



ders. In our case, however, such remainders are simply absorbed into the definition 
of the soft function Sli, through Eq. (^5|). 

The only properties of Sli that we shall need are that it is infrared safe and that 
it has at most a single, infrared logarithm per loop in the limit ws — > 0. Its infrared 
safety follows from the universality of collinear singularities, which have been absorbed 
into the ^'s. Its single-logarithmic dependence on ws requires further explanation, 
however. 

Double logarithms require both collinear and infrared enhancements. Collinear 
logarithms must arise from momentum configurations in which lines in Sli that are 
parallel to the initial-state jets are much more energetic than soft lines that connect 
the two jets. But this requires that the soft lines are much softer in energy than 
wsMjj, since the total, and hence the jet energy is bounded by this quantity. We 



may then apply the reasoning of Ref. |51] to any such momentum configuration to 



show that such a region cancels, until the soft gluons become energetic enough so 
that they may no longer be ignored kinematically compared to lines that are parallel 
to the incoming momenta. But in this region, we cannot generate collinear logs, and 
we expect at most one logarithm per loop. We shall review the explicit forms of the 
soft and jet functions in the following section. First, however, we discuss the role of 
moments. 

3.3 Moments of the partonic cross section 

From Eqs. (|23| ) and (p5|), we find for the full partonic cross section, 
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a s (jjf 



X j(/l) ^ 1; _^ ;as(yU 2 )A j j(/ 2 ) U 2; _^ )as(// 2 ))52 

+0(1 -r), 

where, as above, r = Mjj/ S. In moment space, this becomes 
dvf A f B -+j 1 j 2 (S,5i,5 2 ) _ 

f JZ 



(29) 
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dMjjdAy 



// 



X ^/a//a ^ 



M 



,/,/ 
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a s (/x 2 ),e 



x S 



V /iiV 



A|/,a s (/i 2 
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x JCA) ( iV, c^ 2 ), ) JCA) ( N, a s (^), *. 



M 



j. 1 



+OQ./N) . (30) 

Comparing Eqs. (0) and (0) we derive for the Mellin transform of the hard scattering 
function the "refactorized" expression, accurate to 0(1/N), 

~i>fA/f A (N, Mjj I >, e)#f B/fB (N, Mjj/ /I, e) 



a t {N) 



<t>fA/f A ( N i ^ 2 > e )^/ s // B ( N > ^ 2 > e ) 



7L 





Mjj 
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/i 


(• 






/i 



a s (/x ),5 2 



(31) 



The first factor is "universal" between electroweak and QCD-induced hard processes, 
and was computed first with ]a = <1 for the Drell-Yan cross section || . 

To organize both the \i- and iV-dependences of the jet and soft functions, we 
develop definitions for them in terms of matrix elements in the next section. We will 
then go on to show how the renormalization of these matrix elements leads to the 
resummation of logarithms of N. 



4 Soft and Jet Functions 

We now turn to the explicit forms of the jet and soft functions that result from factor- 
ization, in terms of matrix elements of composite operators. It is the renormalization 
properties of these operators that will lead to resummation of threshold singularities, 
in the next section. 
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4.1 Center of mass distribution 



We begin with the center of mass parton distribution functions ■*/>/// 05 appearing 
in the "refactorized" expression, Eq. (|29f) , where they absorb long-distance contri- 
butions of the initial-state jets, while respecting the overall phase space restrictions 
near partonic threshold. The functions tp differ from standard light-cone parton dis- 
tributions by being defined at fixed energy, rather than light-like momentum fraction. 
They were introduced for the inclusive Drell-Yan cross section, and applied as well to 
heavy quark production |T2" |. For completeness, we may define them by analogy to 
light-cone parton distributions via the matrix elements, 



ip q/g (x,2p /fx,e) 
ip g/g (x,2po/n,e) 



2tt2 3 / 2 



dy e- ix ^(q(p)\q(y , 0) -v ■ 7 q(0)\q(j>)) 



2tt2 3 / 2 
1 

2?r2 3 / 2 p 



dy e" 



-ixpoyo 



(s(p)l Tr 



-v ■ 7 q(y o ,0)q(0) 



k(p)) 



+ , dy e- ixp ^(g(p)\v,F^(y , 0) v v F" ± (0)\g(p)) 



(32) 



where the matrix element is evaluated in n-A = gauge in the partonic cm. frame (in 
Ref. ||, A = gauge was chosen). The vector v is light-like in the opposite direction 
fromp' 1 , so that for p in the ±3 direction, v -7 = 7 ± . In the antiquark distribution, the 
combination qq is treated as a Dirac matrix to define the trace. Charge conjugation 
invariance implies that ip q / q = ipg/q- The moments of these distributions may be 
factorized into a product of moments of the light-cone parton distribution 0, defined 
in any scheme, times an infrared safe function 0. The argument e in tp represents 
the universal collinear singularities that ip absorbs in the factorized cross section, Eq. 



4.2 Matrix elements for final-state jets 

As the second step in the construction of the functions into which the cross section 
factorizes, we treat the final-state jets In axial gauges, all collinear logarithms 

are generated by the imaginary parts of two-point functions. In the factorized cross 
section, Eq. (|29|), the final-state jets are linked to the other functions through the 
weight convolution. The only direct dependence on the momentum of final-state 
particles associated with is in the hard-scattering function Hfl, which depends 
on the total momenta emitted into the cones through Ay and Mjj. In Hjl, the 
final-state (and initial-state) jet momenta are approximated by light-like momenta 
Pi — Pi Mjj/ a/2, with fa a light-like velocity vector in the direction of the jet. The 
vector p^ is determined by the three-momentum of particles flowing into the jet cone 
at partonic threshold. In constructing the jet function, we will therefore sum over 
the complete phase space of the particles associated with jw, subject to fixed total 
spatial momentum of particles within the cone, and fixed weight Wi. 

To be specific, let ii be the total momentum of the final-state particles of 
Pi be the total momentum of particles of J^' emitted into cone i, and ki the total 
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momentum of the particles of J^v emitted outside cone i, so that ii = Pi + fcj. The 
sum over phase space includes the integral over all hi, and over the invariant mass 
p 2 , at fixed values of p^°\ defined as above. The remaining phase space is that of 
two light-like vectors p[°\ one for each of the two jets. These vectors are fixed in our 
differential cross section for Mjj and Ay, Eq. (Iffi) . Overall factors associated with 
this two-particle phase space are implicitly absorbed into the definition of the Hjl- 

Consider, for example, an outgoing quark. Near threshold, the jet function is the 
cut quark two-point Green function, summed over final states at fixed Wi and p^ '. 
This two-point function depends, in general, on p^°\ and on the gauge vector, n M . It 
may thus be defined as 

4«M (P^^ii Mjj, H, a s ^ 2 ), 6 t ) = £ 2^ | (2vr) 3 5 3 fa® - 

C 

X 6(Wi-W^,6i)) (0| / A6 (0) /a,a(0) |0), 

(33) 

where is the field of flavor /, with Dirac and color indices (3 and h, respectively. 
As usual, we suppress dependence on n M . The factor 2\p t (0) \ (2vr) 3 <5 3 -p ? ) in 

Eq. fl3"3| ) matches the normalization of phase space, and fixes p^, the spatial 
momentum of all particles in state £ that flow into the cone. The sum is over all 
states £ with the specified jet momentum, consistent with a contribution Wi to the 
overall weight. 

In our examples, the weight is given by either the first or second line in Eq. 
(p2]). With either definition, as u>j vanishes, the momenta of all particles outside 
the cone are forced to vanish. In the case M 2 j = 2pi ■ p 2 , for example, final-state 
particles from contribute through the term pj/Mjj if they are within the cone, 
and through 2k®/Mjj if they are outside the cone, as in Eq. (|27|). The function 
w {£, $i) — Pi/Mjj + 2k® /Mjj is simply the sum of these two contributions, where 
is the total energy of final state particles in £ that flow outside the cone Si. In our 
approximation, we may neglect corrections to J^ % > that are proportional to Si, the jet 
invariant mass, p 2 , and lOj. 

Near partonic threshold, the invariant mass of jet % is given by 

p 2 i ~2p? ) -p i = V2Mj J f3 i -p i , (34) 

where again fa is a light-like vector in the jet direction. In a frame where (3 is in the 
+3 direction, fa ■ p, is the minus (opposite-moving) light-cone component of the total 
momentum of the particles emitted into the cone. (Recall that by the jet definition 
of Sec. 2.1, p iiT = 0; see Eq. (§.) 

Another consequence of our approximations is that for a quark jet the only Dirac 
matrix structure that we need to retain is proportional to 7 • p^ , 

J ( £!ba {Pi { °\wi, Mjj, fi, a s {fl 2 ), Si) 

= (7 ■ pf ] ) pa ha J {U) U, ¥*L, a s (/i 2 ), S^j , (35) 
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where the function J^*> without indices is a scalar distribution in u>j. This is the 
function in the factorized cross section Eq. (|29|) . Corrections to ( |35l) are suppressed 
by a power of Wi. Similar considerations apply to gluon jets. 



4.3 The eikonal cross section 

Having identified matrix elements that describe initial-state and final-state jets, it 
only remains to give an explicit construction for the soft functions in Eq. (|29|), 
in terms of matrix elements of composite operators. We shall proceed in two steps. 

We start by representing the coupling of soft gluons to the partons involved in the 
hard scattering, of flavors /a, /b, fi and f 2 , by ordered exponentials, also known as 
Wilson lines. We shall introduce the notation 

${/ l) (A 2 ,Ai;x) = Pexp (^-ig drj f3-A {f *\r](3 + x)^ , (36) 

with A™*' the gauge field, represented as a matrix in the representation of flavor 
i = A, B,l or 2, of the gauge group SU(3), and (3 the velocity four- vector of the 
parton whose interactions with soft gluons are being approximated. The operator P 
orders group products in the same sense as the ordering in the integration variable 
A, with the A's with lower values of A to the right. The difference in $^ and $^ 
is in the matrices A® = A%(\ a /2) and A® = A£(-\*J2), with A a the Gell-Mann 



matrices. Note that this notation differs slightly from that in Ref. |12 . 

As we have observed above, these ordered exponentials summarize not only the 
coupling of soft gluons to a single quark or hard gluon line, but also to an entire 
jet connected to the hard scattering by such a parton line fl], |[ |19| . By connecting 
these Wilson lines at a local vertex, we construct an eikonal nonlocal operator, which 
describes the emission of soft radiation, due to both incoming and outgoing hard 
partons. We denote the resulting nonlocal operator as w T , 

4 f) 0% fc} = £ f^- °; s^fo o; x)*,di 

R} 



x 



with the (3i the four-velocities of the Wilson lines that represent the initial- and 
final-state jets. The color tensor (c^p\ describes the couplings of the ordered 

V / d2di,d,Bd,A 

exponentials with each other in color space. 

We use the operator wj to define a dimensionless "eikonal cross section" , describing 
the emission of gluons by the ordered exponentials, 



x(0|T [(4 f) (0)) \ bi} ] \0(!;\T [wf\0) {bi} ] |0) , (38) 

where ^ designates a set of intermediate states, whose contributions to the weight are 
given by w(£, Si). Recalling from Eq. ( p4[) that the final-state jet masses are linear in 
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the momenta of particles emitted into the cones, we may identify the contribution to 
the total weight from state £ as 

«;(£, ft) = ^ (39) 

for Mjj = 2pi -p2, where fc^ is the momentum emitted into the ith jet cone, while k'° 
is the energy emitted outside the cones. Similarly, for Mjj = {jpi +P2) 2 , we have 

2k'° 

MUd = Wj - (40) 

Defined in this fashion, the eikonal cross section contains both collinear and ultraviolet 
divergences, which will be treated by factorization and renormalization below. 

Setting aside the ultraviolet divergences for the moment, we note that factorization 
implies that this cross section is a good picture of the emission of soft radiation as 
a result of the hard scattering, at least for soft quanta that are outside the cones 
of the final-state jets, and not collinear to the incoming partons. Inside the cones, 
or in the directions of the incoming partons, the collinear divergences of the eikonal 
cross section are qualitatively the same as those in the full partonic cross section, 
but quantitatively different in general. The differences are due to the original eikonal 
approximation, necessary to factorize soft emission from the jets. 



4.4 The soft function 

To extract the infrared safe soft function S^} from the eikonal cross section Eq. (|38|) , 
we separate collinear and infrared divergences in the eikonal cross section. Contribu- 
tions from collinear quanta are by construction incorporated in the functions ipi and 
Ji in Eq. (G5f). To include these regions in S LI would be double counting. 



We avoid this double counting by eliminating collinear singularities associated 
with the initial-state and final-state jets from S LI . The observation that makes this 
procedure possible is that the eikonal cross section may be factored into initial-state 
and final-state jets, and a left-over "reduced" soft function, S, in the same manner as 
the full partonic cross section. This, "eikonal", factorization is illustrated in Fig. [|. 
The result is a convolution in the weights of final-state partons associated with each 
of these functions. By analogy to Eqs. ( P7| ) and Q2~8|), these weights are given by 

9b' 9b' 

w' A = ^- v/ B = =^- 

A Mjj B Mjj 

, y/2f3i-k[ + 2k' 1 ° , y/2(3 2 ■ K + 2k' 

Wi = — Wo ~ ' 



Mjj Z Mjj 

Mjj 



9k' 

«4 = -rf- (41) 



for Mjj = 2pi ■ p 2 , and by 



9k' 9V 

w' A = l 4f- w' B = ^ 

Mjj Mjj 
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9k' 9k'° 

1 Mjj 2 Mjj 

nu 

™'s = it 5 -, (42) 

for Mjj = (pi + P2) 2 , where we have used Eq. ([34]) to parameterize the contribution 
of the eikonal jet functions to the corresponding jet invariant masses. The primes 
simply indicate that these variables refer to a convolution for the eikonal, rather than 
the full, cross section. 

The soft function S B j found by factoring exactly the same soft function 

found from the factorization of the full partonic cross section, Eqs. (|2"5|)-(|2"51) above, 



because the soft radiation is insensitive to the internal structure of the jets and the 



hard scattering. Thus we have, by analogy to Eq. (^) and (|29|) 



(f,eik) / wMjj 2 
V fi 



1 

dw' A dw' B dw[dw2dw' s 8 (w — w[ — — w' A — w' B — w' s ) 



c=A,B \ r 1 I d=l,2 \ 



xSS(^,A,,aV)) ■ (43) 

( f ) ( f *\ 

Jin an d iouT are initial-state an d final-state jet eikonal functions, respectively, which 
summarize the dynamics of gluons collinear to the Wilson lines of the eikonal cross 
section. They can be given specific operator definitions. 

Consider first the eikonal distributions for the initial-state jets, i = A, B. The 
phase space for the initial-state eikonal jets is defined by the total energy that they 
emit into the final state (see Eqs. fl4"T| ) and P^)). As a result, their definitions are 
similar to those for the full center-of-mass distributions, Eq. (|32|). 

\ /! J 111 J-00 

x (0| Tr{ f [$Cttt( , -00; y)]T[Q%\0, -00; 0)] } |0) , (44) 



with y u = (y , 0) a vector at the spatial origin. As in Eq. (P2F), e in the arguments of 

(f) 

Jin denotes collinear singularities. 

Similarly, the collinear dynamics of the eikonal final-state jets are summarized by 
matrix elements analogous to Eq. (|33| ). 



x <0| Tr { f[^(oc, 0; 0)] |£) (£|T[*J« (00, 0; 0)] j |0) , (45) 
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with i — 1, 2. Here again w(£, Si) is given by either Eq. ( |4l| ) or 

We construct (moments of) the soft function by dividing the moments of the 
eikonal cross section ( |43| ) by the product of moments of the eikonal jets, Eqs. 
and (HD, 



JIN 1, Nfj, ■ 



(46) 



This simply corresponds to the standard factorization procedure in axial gauge [|IJ . At 
one loop, all diagrams that are two-particle reducible by cutting incoming or outgoing 
eikonal lines with the same four- velocity $ in the amplitude and its complex conjugate 
are eliminated from the soft function, which is then free of collinear divergences. 

Both the eikonal distributions Eq. ([44D and (f45"D , and the soft function Eq. (|46|) , 
are as yet unrenormalized. Eikonal lines, and vertices made from their products, 

This 



may be renormalized in the usual multipicative manner f22| , [23| , [24 
renormalization will play an important role in resummation below. 
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5 Resummation 

We are now ready to use the renormalization properties of the jet and soft functions 
to organize iV-dependence in the factorized cross section in moment space, Eq. (|3l"l) . 

5.1 Renormalization of the soft function 

The soft function S^j, as emphasized above, requires renormalization as a composite 
operator. This renormalization is a direct consequence of the factorization we have 
discussed. The ultraviolet divergences of S L j arise when approximations appropriate 
for soft gluons are extended to all momenta. The renormalization scale for the soft 
function, then, acts as an effective cutoff, separating soft from hard gluons. Indeed, 
in the product H IL Su, ultraviolet divergences induced by the factorization cancel 
against each other by construction, since the original diagrams have no UV diver- 
gences beyond those taken into account by the usual renormalization of the theory 

Because H and S occur in a product, they must renormalize multipicatively, with 
separate renormalization factors for the amplitude and the complex conjugate |23|, [7| 



H^fl = jljT l {zf^) io Hc D [{zf^\ DL 
5 (f)(o) = \zf) LB S BAZ V, (47) 
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where Zj is the renormalization constant of the ith incoming parton (flavor ■ ■ ■ f-i) 
connecting to and 

Zscd i s a m atrix of renormalization constants, describing 
the renormalization of the soft function. 

From Eq. fl4"T|), the soft function S^j satisfies the renormalization group equation 

|2T 



s 



(f) 

LI 



(4 S )LBO BI 



(f) 



b LA\ l S JAI 



(4£ 



(f) 

where we have introduced a soft anomalous dimension matrix, Ig , which is computed 
directly from the UV divergences of the soft function. We can compute the matrix in 
a minimal subtraction renormalization scheme, taking e = ejjv = 4 — D, with D the 
number of space-time dimensions. One-loop anomalous dimensions are then given by 



r ( s ) (9) = - g ^Res^ zP(g,e) 



(£) 

The determination and analysis of T s has been carried out in the case of heavy quark 
production in [[EJ; its calculation for massless quarks |23[] and gluons at one loop will 
be the subject of a companion paper. In any case, the solution to Eq. fl48|) takes the 
form 



(49) 
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I V V- 
xPexp[/ AW ^r«t( as(/i 

xS®(l,Ay,a s (M^/N 

M Jj/N dfi' 

' S 



xPexp 



12- 



(50) 



where the trace is taken in the space of color structures. The symbols P and P 
refer to path-ordering in the same sense as the integration variable fi' and against it, 
respectively (for example, P orders Tg\a s (Li 2 )) to the far right and T^p (a s (Mjj /N 2 )) 
to the far left). As usual, at leading logarithmic accuracy in /i, we can simplify this 
result to a sum of exponentials, by choosing a basis in which the matrix T s is 
diagonal. 



5.2 Initial-state jets 

Starting from Eq. ([H]), we recall the resummed expression for the prefactor {ip/<p) 
in Eq. ([31]), the exponentiation of whose iV-dependence follows, as usual, from its 
factorization properties |3], [/j. The general expression for moments of the ratios of 
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the functions ip and 0, evaluated at the common scale /i = Mjj, is 



I 
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(51) 



where 
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(52) 



The parameter ms and the resummed coefficients depend on the factorization 
scheme, that is, on the definition of 0///- The results are rather different for DIS and 
MS schemes, in particular {|, |?J. This difference must be compensated for by differ- 
ences in the parton distributions themselves. With DIS (MS) factorization schemes, 
we have m s = 1 (0) in Eq. (|52|). In Eq. (^l|), R(f)(a s ) is an iV-independent function 
of the coupling, which can be normalized to unity at zeroth order. 

The , B^f) and are finite functions of their arguments. To reach next-to- 
leading order accuracy in In (AT), we need 
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A (f) (a s ) = C f i^ + ^K -± 



7T 



7T 



(53) 



with Cf = Cp {Ca) for an incoming quark (gluon), and with given by |2S 

'67 7T 2 \ 5 



lc c 



9 



n f 



(54) 



where n/ is the number of quark flavors. Bin is given for quarks in the DIS scheme 
by 

B (q) (a s ) = -~C F ^ , (55) 



7T 



while it vanishes in the MS scheme for quarks and gluons. (The DIS scheme is 
normally only applied to quarks, although extended definitions for gluons are possible 
|29fl.) Finally, the lowest-order approximation to which is scheme- independent, 



is 



,(/) 



(56) 



The results just discussed are useful, but not quite adequate for our needs, because 
they assume that // = Mjj, as would be appropriate for Drell-Yan cross sections. For 
jet cross sections, we generally need additional freedom to choose the factorization 
scale. To change the scale /i, we need the renormalization group behavior of the 
parton distributions ip and 0, whose ratio enters Eq. (]3"T|). 

The center of mass distribution ip requires no overall renormalization as a com- 
posite operator ||. From its definition, Eq. fl3~2f ), we see that ip, and each of its 
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moments, renormalizes multiplicatively, because it is the matrix element of a product 
of renormalized operators. As a result, it obeys the renormalization group equation 



d^ f/f (N,Mjj/^e) 
d\i 



2 7/ (a s (/x 2 ))^ /// (iV,M JJ //i,e) 



(57) 



with 7j the anomalous dimension of the field of flavor /, which is, of course, indepen- 
dent of N. 

The dependence of the light-cone distribution on the factorization scale /i is 
slightly more complex than for ip, and depends on the factorization scheme that we 
choose. The simplest case is the MS scheme. Each of the moments of the MS version 
of cf) obeys a renormalization group equation with the anomalous dimension of the 
color-diagonal splitting function for that flavor, 



// 



d4>f/f(N, n 2 
d\i 



2 7// (iV,« s ( / i 2 ))0 /// (iV,/i 2 



(5? 



For the MS distribution, this relation holds by definition, because the distributions 
are defined as the matrix elements of operators on the light-cone, whose renormal- 



ization is described by the splitting functions |3(J, 131 ■ O n ly the flavor-diagonal 
evolution survives in the large- N limit, because only color diagonal splitting functions 
are singular for x — *■ 1. 

The factorization scale dependence now may be controlled, by the two evolution 
equations (|57]) and (|5"8D . The MS scheme expression that includes the /^-dependent 
prefactor in Eq. fl3"T|), thus generalizing Eq. (|5T]), is 
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R 



(/) 



(q s (]U 5 



exp 



E (f) (N,Mjj) 



J MS 



x exp 



M, 



d/i' 



7/ (a s (^))- 7 //(iV,a s (/i /2 )) 



(59) 



This expression enables us to change factorization scales in the resummed cross sec- 
tion. The corresponding factor for DIS scheme distributions is easily found by com- 
paring the ratio in Eq. (|5lD in DIS and MS schemes. The effect of the DIS scheme is 
to produce an extra factor in Eq. (|59f) , 
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<Pf/f( N ^ 2 ,£) J DIS 
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x exp 
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(60) 
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This difference may also be formulated in terms of a slightly modified evolution equa- 
tion for DIS distributions, which takes into account their logarithmic moment depen- 
dence 0. 



5.3 Final-state jets 

The remaining factors in Eq. ( |3~ll) are the final-state jet functions J^ l \ The resum- 
mation of their iV-dependence depends critically on the definition of the dijet cross 
sections. For illustrative purposes, we have identified the choices Mjj = (pi + p 2 ) 2 
and 2p 1 ■ p 2 . We begin with the latter. 

The case Mjj = 2pi ■ p 2 . In this case, the jet function contributes to the overall 
weight through the mass of the cone-jet, p 2 , and also through the total energy of 
particles emitted outside the cone. Thus, recalling Eq. fl27|) , we may write as an 
integral over these variables at fixed Wi, 



—^,a s {n),5i\ = Jdpid^dlwi — 2 

2 uo \ 



W^A*I, (61) 

\ fl /I 2 /I / 



with 1^' a density in the invariant mass inside the cone, and energy outside. This 
makes its analysis a bit more complicated than the cases studied in Ref. 0. 

In the limit of very small pf, however, we may factorize soft gluons emitted at 



angles much larger than yp 2 /Mjj from a function that describes the dynamics of 
collinear partons. This factorization can be carried out by the procedure of Ref. ||, 
referred to in Sec. 3.2 above. As discussed above, couplings of soft gluons to the jet are 
first replaced by an eikonal, or "soft" approximation. Suppose the jet velocity is 
In the soft approximation, soft gluon momenta carried by jet lines are approximated 
by their "opposite-moving" component Pi-q, and similarly for soft gluon polarizations 
at the vertices where they are emitted by jet lines. 

After the use of Ward identities, soft gluons factorize in the soft approximation 

H| , and we derive a convolution of the form, 



Mjj p] fc? . \ I JIO , (Mjj v \ - V2Mjj(P, ■ q) 



\ u u z u J 



/l '11* pi J J \ jJ, /i 2 

;(fi) 




^F^.- • (62) 



Corrections to this relation are free of infrared divergences in the limit of small kf, 
(Pi ■ q), and hence do not contribute to the leading power for tUj — > 0. The function 
iouT' which absorbs the soft gluon dynamics is the same eikonal jet function defined 
in Eq. (|4"5| ) and shown in Fig. |J The function remains dependent on the jet's 
longitudinal momentum, and is linked to Jqut by a convolution in the small compo- 
nent (/3i ■ q) of momentum emitted by that function into the cone. In the limit of 
small p 2 , gluons emitted outside the cone are all associated with the eikonal function. 
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The most important feature of S'-W for us is that it is independent of the cone size, 
up to corrections that vanish with p 2 . 

The convolution in Eq. (|62| ) simplifies under moments of Eq. (pT]), with 

respect to u>j, 



V A* 



N { ^TT- J.c/0 / (A • g) M 



/i /i y 

where in the second line, p' 2 = p 2 — \f2Mjj((3i ■ q). The moment integrals here and 
below all have lower limits at zero. Dependence on their upper limits is exponentially 
suppressed in N, and may be neglected. Our task now is to derive the iV-dependence 
of the product j^t^l 

To find the requisite iV-dependence, we begin by observing that a very similar 

procedure may be applied to the full axial gauge two-point function, which we denote 

( f ) 

as £ 2 • Again, soft gluon emission away from the jet axis factorizes into a convolution 
form, 



pi ' J J V A 2 ' ^ 2 J V A 4 

'(64) 

f f) 

with j_ the two-point eikonal function, with its eikonal line again in the /^-direction, 
evaluated at fixed values of the /?$ • q momentum component of emitted partons. Here 
is the same function as in Eq. (|62l) for small p 2 , because in this limit collinear 
gluons are emitted only at the center of the cone, far from its boundary. Just as for 
Eq. ( |63|) , £2 factorizes into a product under moments, this time with respect to pf, 



x 



Comparing Eqs. ( |63"D and (p5]), we see that the moments of the cone-jet function Jv*) 
with respect to Wi are closely related to the moments of the full two-point function 
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S 2 with respect to pf, 
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(66) 



The A/"- dependence of the moments of the jet function are thus given by the moments 
of the full two-point function, times the ratio of moments of the two eikonal functions 
Jout and j_. We readily find their iV- dependence as follows. 

Moments of a cut two-point function with respect to its invariant mass were re- 
summed explicitly in Ref. ||. The ln(iV) dependence exponentiates in analogy to 
Eqs. ( |51~1) and (|52|) . Normalizing the two-point function to 5(pf/Mjj) at lowest order, 
we have 
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where 
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(68) 



The prime on E' indicates the exponent for a final-state jet, while the function Am 



The lowest order term in B',f\ may be read off from the 



is the same as in Eq. (| 
one-loop jet function. The results include a gauge dependence, which cancels against 
a corresponding dependence in the soft anomalous dimension matrix [23]. 

The other factor in Eq. (|66"D , Jout(A0 / actually affects the moments of 
the final-state jet only beyond next-to-leading logarithm in N. This may be seen as 
follows. First, being the sum of fully eikonal diagrams, both j'out and j- exponentiate 
in moment space, according to the general arguments of Ref. [32|. At the same time, 
double logarithmic contributions, which are associated with gluons collinear to the 
eikonal lines, match in j'out and j_, because the phase space for collinear gluons 
is the same in both functions. The ratio may thus include at most next-to-leading 
logarithms, associated with soft gluons emitted outside the cone. The contribution 
of next-to-leading logs in the exponential may then be determined from a one-loop 
calculation. As for the soft function above, single logarithms are determined by the 
renormalization of one-loop virtual diagrams, which are identical for j'out and j_. 
An explicit calculation therefore gives no logarithms at all in the ratio. 

In summary, up next-to-leading logarithm, Jout/J- — 1> and the moments of the 
final-state jet in Eq. (|66|) with Mjj = 2p 1 ■ p 2 are given by 



d Wi e- Nwi J {h \wi) 



exp 



jj, 



(69) 



with E',f\ defined as in Eq. (|68D, up to next-to-next-to- leading corrections in the 
function B',^. 
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An important feature of the exponent associated with the final-state jet is that 
it has the opposite overall sign from the exponent of Eq. ([52]) , associated with the 
initial state. The leading logarithmic contributions of the final-state jets are negative, 
but are otherwise determined by the same functions A^(a s ). Just as the leading-log 
initial-state contributions always act to enhance the cross section, those associated 
with the final state always suppress it, by an amount that depends on the partonic 
subprocess. Such singularities correspond to Sudakov suppression of scattering in 
the elastic limit |7], 2~3fl . They are already present in explicit next-to- leading order 



calculations of single-jet inclusive cross sections [331 - It ma y be worth pointing out 
that for the specific choice of subprocess qq — > gg, the net coefficient of the leading 
logarithm is negative in both DIS and MS schemes (because Ca > Cp/2), in sharp 
contrast to the classic Drell-Yan calculation |3|. [| . 

The case Mjj = (p\ + P2) 2 ■ When we choose Mjj as the total invariant mass 
of the two-jet system, the phase space at partonic threshold is, as mentioned in Sec. 
3.1 above, changed significantly compared to the previous choice. In the present 
case, the masses of the dijets are not forced to zero at threshold (see Eq. fl2"2|)). 
In fact, the phase space at partonic threshold consists of all states in which each 
of the two jets carries total energy Mjj/2, with equal and opposite momenta. By 
dimensional considerations, the jet masses may increase to the order of 5 2 Mjj. As a 
result of this larger phase space, collinear logarithms of N cancel, leaving only infrared 
enhancements associated with soft emission outside the cone, which appear as single 
logarithms of N per loop. The remnants of the dilogarithmic structure near threshold 
are terms of the form a™ (ln<5; In N) n , which replace the leading a™ (In N) 2n terms 
of the previous case. 

Although the jet functions with Mjj = (pi + p 2 ) 2 involve fewer logs per loop 
than for Mjj = 2pi -p 2 , the situation is in some sense more complicated than before. 
This is because of configurations in which the jets, although of minimum energy, are 
nevertheless massive. In these regions of phase space, soft emissions near the edge 
of the cone may resolve the flow of color within the cone []14 , carried, for example, 



by fast partons separated by angles of order Si. The couplings of soft gluons to such 
configurations may not be approximated by their coupling to a single eikonal line at 
the center of the cone. This approximation was accurate to all logs for Mjj = 2pi-p 2 , 
because in that case the limit N — > 00 forces the jet mass to zero, and hence the 
energetic partons to the center of the cone. Nevertheless, the approximation in which 
soft gluons are emitted by a single eikonal line at the center of the jet still captures 
the next-to-leading, a™ In™ N, contributions. Any such eikonal cross section, with 
phase space that is symmetric in all particles, exponentiates p2| , [34], 



(V, a s (/rV*l = exp[E' (f) (N,Mjj)} 



E{ f) (N,Mjj) = exp 



J (f) 

/ v c u)Mv )) 



(70) 



where C'^(a s ) is a finite perturbative series whose first term may be read off from a 
one-loop calculation. Beyond these next-to-leading logarithms, however, corrections 
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to Eq. ( |T0D give a rather complicated sum of terms, each representing the coherent 
radiation due to a set of eikonal lines along arbitrary directions within the cone. 

Again to avoid a proliferation of symbols, we use the same notation for the expo- 
nent here as for Mjj = 2p\ -p2, although the leading behavior in the two cases is quite 
different. For Mjj = (p% + P2) 2 , the leading logarithmic behavior in moments (and 
therefore the leading distributions in momentum space) are not affected by the final- 
state jets, and thus retain the same (positive) contributions encountered in Drell-Yan 
cross sections. 

We are now finally ready to assemble all the pieces necessary to write down a 
resummed cross section for dijet production. 



5.4 The resummed cone-dijet cross sections 

The characteristically nonabelian aspect of resummation in our case lies in the evo- 
lution of the soft functions. Recall that because the anomalous dimension matrix is 
not diagonal in general, solutions to the evolution equation ( |i8|) are ordered, rather 
than simple, exponentials. Substituting the solution (|50f) to Eq. (f|8|) into the cross 
section in moment space, Eq. (pi]), and using Eqs. (|5TD, (|59D and (^) (or Eq. (|70|)), 
we find an expression that organizes all logarithms of N at leading power in N, in 
MS factorization scheme, 
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(71) 



where as above the trace is taken in the space of color structures. For DIS scheme, 
the only change is the additional factor given in Eq. (]60f) . 

Without going into detail at this point, certain general features of the resummed 
cross section may be readily identified. The first exponential factors, exp[Erf.\(N)], 
serve, as in Drell-Yan cross sections, to enhance the cross section. The next factors, 
involving the anomalous dimensions of the parton fields and splitting functions, gov- 
ern factorization scale dependence. The third factor, associated with final-state jets, 
always acts to suppress the cross section. We have noted already how the choice of 
jet algorithm can influence their size. The hard-scattering functions remain depen- 
dent on the scattering angle. The most interesting factors, however, are the ordered 
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exponentials of the soft anomalous dimension matrix. They depend on the scattering 
angle and they distinguish the roles of different color structures in the hard scattering 
through coherent gluon emission, linking initial- and final-state partons. 

The logarithmic accuracy of the resummation depends upon the order to which 
the hard scatterings, and the anomalous dimensions, have been calculated. In a paper 
to follow this one, we shall compute the leading-order anomalous dimension matrices 
for each flavor combination in QCD. This will enable us to present explicit resummed 
dijet cross sections to next-to-leading order in ln(iV) [ |35|j . 



5.5 Other jet algorithms and kinematics 

Eq. (|TT|) is typical of resummed jet cross sections, but its details depend on our use of 
cones to define the jets. Other choices, of course, are possible, and in cases, preferable. 
One of the disadvantages of the cone criterion is that it requires an extra parameter 
Si for each jet. In addition, single-jet, rather than dijet cross sections are often more 
convenient experimentally. Single-jet inclusive cross sections have slightly different 
kinematic properties near threshold, compared to dijet, and therefore result in slightly 
different resummed expressions [13 . The extension of threshold resummation to cross 



sections with single-particle kinematics will be discussed in Ref. [[36]. Nevertheless, 
it is clear that the result ([H]) is much more general than the cone algorithm that 
we have used. In fact, all that is necessary to derive such a resummed expression 
is a convolution form like Eq. (^), in which the weight is linear in the momenta of 
radiated gluons. Threshold logarithms in any such weight will be controlled by the 
same matrix of anomalous dimensions T$ that we have identified above. Differences 
between different algorithms will, in general, show up only in the functions B',^ for 
the final-state jets in Eq. (|68]). 



6 Conclusion 

In this paper, we have shown how to resum threshold logarithms in cone-based dijet 
cross sections. The method may be extended to other jet algorithms, and the explicit 
form of the resummed cross section will depend upon the details of the cross section 
chosen. In every case, however, the resummation depends upon the color structure of 
the hard scattering. As for the Drell-Yan cross section, the resummation of threshold 
singularities depends on the factorization scheme, although the color- dependent con- 
tribution does not. Our explicit result for the MS factorization scheme is given in Eq. 
(fn]) above. It is, in principle, valid to all logarithmic accuracy, at leading power of N 
in moment space for the dijet cross section at fixed p\ -p2, and to next-to-leading log- 
arithm in N for fixed (pi +P2) 2 - Its inverse transform to momentum space therefore 
summarizes singular distributions in 1 — z to all logarithmic order in the former case 
and to next-to-leading order in the latter. As for heavy quark production, the general 
resummation can only be given in terms of ordered exponentials, due to mixing of 



color exchange tensors by soft gluon emission. In a forthcoming paper [35], we will 



derive the soft anomalous dimension matrix to one loop for the full range of flavor 
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scatterings that give rise to jet production, and give explicit expressions for dijet 
cross sections to leading order in the coupling, and to leading and next-to-leading 
logarithm in the moment variable. 
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/P 2 

Figure 1: Configuration of cones in the partonic center of mass frame at partonic 
threshold. 
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Figure 3: Factorized form of the cross section into initial-state and final-state jet 
functions (ip and J, respectively) and the color-dependent soft function S^j . At lowest 
order, S L j is given simply by the set of all diagrams in which a single gluon connects 
any two ordered eikonal lines moving in different directions. The construction of S L j 
beyond lowest order in discussed in Sec. 3.2. 



34 




Figure 4: Factorized eikonal cross section. With an appropriate choice of color nor- 
malization factor I/A/", the soft function S L j is the same as in Fig. |3|. The eikonal jet 
functions Jout and jw are defined in Eqs. (f|3|) and (fJ5|), respectively. 
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